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ABSTRACT 

Cold gas entering the central 1 to 10^ pc of a galaxy fragments and condenses into clouds. The 
stability of the clouds determines whether they will be turned into stars or can be delivered to the 
central supermassive black hole (SMBH) to turn on an active galactic nucleus (AGN). The conventional 
criteria to assess the stability of these clouds, such as the Jeans criterion and Roche (or tidal) limit, 
are insufficient here, because they assume the dominance of self-gravity in binding a cloud, and neglect 
external agents, such as pressure and tidal forces, which are common in galactic nuclei. We formulate 
a new scheme for judging this stability. We first revisit the conventional Virial theorem, taking into 
account an external pressure, to identify the correct range of masses that lead to stable clouds. We 
then extend the theorem to include an external tidal field, crucial for the stability in the region of 
interest - in dense star clusters, around SMBHs. We apply our extended Virial theorem to find the 
correct solutions to practical problems that until now were controversial, namely, the stability of the 
gas clumps in AGN tori, the circum-nuclear disk in the Galactic Center, and the central molecular 
zone of the Milky Way. The masses we derive for these structures are orders of magnitude smaller 
than the commonly-used Virial masses (equivalent to the Jeans mass). Moreover, we prove that these 
clumps are stable, contrary to what one would naively deduce from the Roche (tidal) limit. 

Suhjeet headings: hydrodynamics - instabilities - ISM: kinematics and dynamics - galaxies: active - 
galaxies: nuclei 


1. INTRODUCTION 

An essential ingredient in galaxy formation and evolu¬ 
tion is the impact of the growth of superm assive black 
holes (SMB Hs) on the galactic environment (Kormendy 
& Ho |20I3 ). The agent is gas - a small amount of cool 
gas with low angular momentum fuels the central SMBH, 
turning on an active galactic nucleus (AGN), which by 
means of radiation and outflow heats the gas content on 
larger, galactic scales, hence regulating the star forma¬ 
tion rate in the entire gala xy and possibly cutting off 
the gas supp ly to the AGN (iMcNarnara & Nulsen||2007 
Fabian|[2^ . 


A bottleneck occurs when cold gas is funneled to the 
nuclear region of a galaxy and forms a disk (because of 
angular-momentum conservation), at a distance of a few 
hundreds to several parsec (pc) from the SMBH. At this 
stage, the gas density rises to the critical point of trig¬ 
gering the Jeans instability, and hence the gaseous disk 


fragments and becomes clumpy (Shlosman et al. 1990, 
Goodman|2003 ). From this point on, the inflow of gas is 


carried on mainly by clouds. But it is known that a cloud 
produced by the Jeans instability is self-gravitating, and 
will collapse on the free-fall timescale. The consequence 
is that within a short time most clouds would turn into 
stars, so essentially, and contradictory to man y observa¬ 
tions, little gas is left to fuel the central AGN (Goodman 


& Tan 2004 Tan & Blackman 2005 Collin & Zahn|20U8 T 
lb resolve this contradiction, several ideas have been 
proposed. In particular, two possibilities to deliver gas 
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clouds to the AGN in spite of them being Jeans unsta¬ 
ble have been considered, (i) A cloud can b e stabilized 


by a strong magnetic fie ld or fast rotation (Rees 1987 
Vollmer &: L)usclil|2001b ). (ii) The inflow of gas could be 


accelerated, as a result of an asymmetric gravitational 


potential 1 

Shlosman et al. 11990 

Macieiewski et al.|2002|) 

or cloud-c 

ioud collisions (IKroll 

k & Begelman 

1988 Ku- 

mar|1999 

Vollmer et al.|2004;|). 

so that gas clouds can be 


stars. 

An alternative viewpoint is that the gas clouds, after 
the initial fragmentation, have gone through a series of 
physical changes, so that the eventual physical proper¬ 
ties no longer satisfy the Jeans criterion. Two physical 
processes could cause such changes: (i) Big clou ds could 


break into s maller ones due to rapid cooling (Burkert 
& Lin |2000). Each descendant cloud can be stable be- 


cause of the reduced self-gravity. A galactic nucleus is an 
ideal place of producing such clouds, because the high 
gas density there enhances the cooling rate, (ii) Colli¬ 
sions between gas clouds are frequent in galactic nuclei, 
since the number density of cloud is high an d the orbital 
periods are short (Krolik & Begelman||l988 ). Such colli¬ 
sions raise the internal pressure of a cloud, supporting it 
against gravitational collapse. 

It is possible today to test the stability of gas clouds 
with observational data, (i) Clumps have been seen in 


and active I 

dunt & Malkan 2004 

Mazzalay et al. 

2013 

Davies et a 

.2014). (il) I'or AG!J 

s, it has been widely 


eral pc must exist surrou nding the SMBH to cause the 
type I/type I I dichotomy (Krolik fc Begelman|1988 An- 


tonucci 


19931. Indeed, direct evidence of such clumps has 


been found by modeling the X-ray variability of AGNs, 
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so that we have a way of measuring some of the clouds’ 
physi cal properties (Risaliti et al. 2002 Markowitz et al. 
2014 and references therein), (iii) In the M ilky Way, 
about 200 molecular clumps have been found (Miyazaki 
& Tsuboif 2000 Oka et al. 2001) within a distance of' 
200 pc from the Galactic Center (GC), in th e region 
known as t he “central molecular zone” (CMZ, Mezger 


et al. 


1996). Further inside, at a distance of 1 to 3 pc 
from &gr A* (the SMBH, 8 kpc from Earth), there is the 
circum-nuclear disk (CND), which has been resolved into 


cular 

dGenzel et al. 11985 IJackson et al.|1993 

Marr et al. 

1993 

Marshall et al. 

TU95). 


terion is required. However, the conventional Jeans cri¬ 
terion and Roche limit become insufficient, because they 
regard self-gravity as the only force binding a gas cloud, 
which is not true in galactic nuclei. Here gas clouds 
are likely subject to a compression by a high external 
pressure. This pressure exists because the interstellar 
medium (ISM) in the nuclear region of a galaxy usually 
has a high density and a large turbulent velocity. More¬ 
over, there is likely a strong tidal force acting onto the 
clouds, because of the ubiquity of SMBHs and dense nu¬ 
clear star clusters (see Amaro-Seoane|2(n2 for a review). 


Here we show that these agents are crucial in determin¬ 
ing whether a cloud is stable, as well as in the correct 
derivation of its mass. 

The paper is organized as follows. In Section we 
revisit the Virial theorem so as to link the physical prop¬ 
erties of a gas cloud to its environme ntal conditions, in¬ 
cluding the external press ure ( Section [2.1[ ) and the back¬ 
ground tidal field (Section |2. 2 [ ). After taking thes e ex ter- 
nal factors into account, we formulate in Section [2.3| our 
''''extended Virial theorem”, which gives only one stable 
solution for the mass, and we discuss its implications. In 
the light of this stable solution, we proceed in Section^to 
interpret the observational data of gas clumps, incl udin g 
those in AGN tori (Se ction [3T] ) , the GND (Section |3.2[ ), 
and the CMZ (Section |3.3[ ) of the Milky Way. Finally in 
Sectionwe justify the assumptions that we adopted in 
the wor^ and we discuss the implications of the results. 


where p is the fluid density, r is its position vector, the 
dots denote time derivatives, P is the internal pressure 
including both thermal and turbulent motion, and is 
the gravitational potential (self-gravity). If we take dot 
product with r on both sides of the last equation and 
integrate it over the entire volume of the cloud, the left- 
hand-side (LHS) of Equation ([^ becomes 

J^p{r-r)dV=^^J^pr^dV-J^v^pdV. ( 2 ) 

If the internal structure of the cloud is steady, the first 
term on the right-hand-side (RHS) of Equation ([^ van¬ 
ishes because is constant. If the cloud is also not 
rotating, the second term vanishes, because v, the macro¬ 
scopic velocity of a fluid element, is zero. 

Since we are integrating Equation (IH, when the inte¬ 
gral of the LHS is zero, the integrated RHS must also be 
zero: 


0 = - / (r • VP -I- pr • V$c) dV (3) 

Jv 

= [ [PV • r - V • (rP)] dV- [ {r- V$c)dm (4) 
Jv Jv 


= / 3PdV — / Pgr ■dS+ r • SL^^-dm, 


( 5 ) 


fV 


where Ps denotes the pressure at the surface of the 
cloud, dm = pdV is the mass of the fluid element, dS 
is the differential area vector orthogonal to the surface 
of the cloud, and we have introduced for convenience 
agi- = — V$c, the self-gravity of the cloud. We now as¬ 
sume spherical symmetry and homogeneity for the cloud 
to simplify the integration. We note however that this as¬ 
sumption does not have an impact in the general conclu¬ 
sions that we will draw. Homogeneity leads to a constant 
density p and a constant one-dimensional velocity disper¬ 
sion a inside the cloud, so P = pcr^ is also constant. In 
addition with the assumption of sphericity. Equation ([^ 
reduces to 


2. VIRIAL THEOREM 
2.1. External Pressure 

A strong external pressure affects the motion of the 
fluid at the surface of a gas cloud, and changes the condi¬ 
tion of stability of the entire cloud. Such a high pressure 
is found in galactic nuclei because of a concentration of 
hot ionized- an d turbulent cold gas. The Virial theorem 
(Clausius 1870) accountin g for an external pressur e has 
been derived in the past (Spitzer||1978| |Shu||1992|) and 
has been applied to the molecular mouas m the (JIVIZ of 
the Milky Way (e.g. Miyazaki & Tsuboi||2000 Oka et al. 
2001). In this subsection we present the derivation m 
detail so as to highlight an important list of conceptual 
points - which are often overlooked in the literature - 
that will be crucial for the main ideas of this work. 

Without magnetic fields and tidal forces affecting it, 
the equation of motion of a fluid element can be written 
as: 


AttRIPs = 3Ma^ - aGM^/Rc 


( 6 ) 


pr = —VP — pV$c 


( 1 ) 


the conventional Virial theorem for a steady, non¬ 
rotating cloud embedded in a pressurized medium. Here 
Rc is the radius of the cloud, M is the cloud mass, G is 
the gravitational constant, and a is a geometrical factor 
of order unity, which is 3/5 assuming sphericity. 

The connection between the above Virial theorem and 
the external pressure Pext comes through the surface 
pressure Pg - a cloud in hydrostatic equilibrium with 
the surrounding medium will have Ps = Pext- When 
Pext = 0, Equation has only one non-trivial solu¬ 
tion Mg := M = 3Ra^{aG), that is the commo nly-used 
virial mass in the literature (e.g., Zwicky 1937). It can 
be easily verified that Mq is essentially the same as the 
Jeans mass Mj = Xjp (differ by only a factor of 1.02), 
where Xj = a^J-njifip) is the Jeans length. 

In a more general situation, which is what defines the 
fundamentals of this research, Pext > 0, and then two 
non-trivial solutions exist for M in Equation (O, as is 
illustrated in Figure The two solutions M± (and we 
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Fig. 1.— Surface pressure Pg as a function of mass M for a spher¬ 
ical, virialized cloud with constant radius Rc and internal velocity 
dispersion a. Given an arbitrary external pressure Pext between 
0 and the characteristic pressure Mocr^/(4y), two solutions ex¬ 
ist for M that allow the cloud to be in hydrostatic equilibrium, 
i.e. Ps = Pexp- One solution, M+, which is unstable (see text), 
falls in the mass range between the commonly used virial mass, 
Mq = 3Ro-^/{aG), and half of it, Mol2. The other solution, M_, 
is stable and has a maximum value of Mol2. 

note that we adopt the notation of |Oka et al.||200H with 
M_|_ > M_), can be expressed in a compact format 

M± = a±Mo, with a± := - ^ (7) 

where we have defined the dimensionless mass a := 
M/Mq and the dimensionless surface pressure (3 := 
APsV /where V = 47ri?^/3 is the volume of the 
cloud. 

Thanks to these dimensionless parameters, we can 
rewrite Equation (§ in a dimensionless way: /3 = 4a — 
4a^, which is useful to understand the role of these pa¬ 
rameters in assessing the correct values for the mass of 
the cloud, as we will see later. Figure displays the 
geometric meaning of these dimensionless parameters. 

A close analysis of Equation (1^ reveals crucial infor¬ 
mation about the stability of gas oouds in galactic nuclei: 

1. Given (for fixed) i?c and tr - quantities that can 
be extracted from observations - the maximum Pg 
that a cloud can have is Mo(t^/( 4V). This is also 
the maximum external pressure that a cloud can 
withstand in hydrostatic equilibrium. At this crit¬ 
ical point (because it is an inflection point, see 
Figure [^, Pext = Moa'^/{4V), the only solution 
to maintain hydrostatic equilibrium of the cloud is 
M = Mo/2. 

2. When 0 < Pext < Afo(T^/(4V), we have that 
Mq/2 < M+ < Mq and 0 < M_ < MqI2. Hence, 
in this regime of Pext, Afg is always an overestima¬ 
tion of the real mass of the cloud. 

3. Only one branch of the solution for the mass of 
the cloud is dynamically stable. To see this, let 
us consider a cloud already in hydrostatic equilib¬ 
rium, and then we perturb it by either increasing 
or decreasing its mass, (i) First branch: M = M+. 



Fig. 2. — Same as Figure^but plotting the dimensionless quan¬ 
tities. The mass M of a cmud is normalized by the commonly- 
used virial mass Mq, so that cx. = M/Mq. The surface pressure 
Ps is normalized by the characteristic pressure Mocr^/(4V), so 
d = 4PsV/(Moa^). 

If we slightly increase M, the surface gravity of 
the cloud increases, meanwhile the surface pressure 
of the cloud decreases because dPg/dM < 0, as 
shown in Figure The tendency of the cloud is 
to shrink, but this will increase the surface gravity 
even more, so the cloud collapses. On the other 
hand, if we slightly decrease M, then the surface 
gravity decreases while Pg increases, so the cloud 
will expand, which leads to an even smaller surface 
gravity, therefore the cloud eventually explodes. 

(ii) However, for the second branch of the solution, 
M = M_, we have dPg/dM > 0, so the cloud is in 
the opposite situation - Pg changes in proportion 
with the change of the surface gravity. Now the 
cloud is dynamically stable. 

4. As a consequence of the last point, we can see that 
only M_ is physically correct. Observationally, gas 
clouds are often found in places where the observed 
pressure (i.e. not derived from the Virial theorem) 
is such that Pext ^ Moa'^ j{4V), or ,5 1 in the 

dimensionless expression, as we will see later in the 
practical examples. Hence, a common practice in 
the literature is to use Pg = 0 as an approxima¬ 
tion (without any justification a priori). As a re¬ 
sult, the only non-trivial solution of Equation (1^ 
is M = Mq, which is the commonly-used virial 
mass. Nonetheless, we know that there are, strictly 
speaking, two solutions for M even when Pg is very 
small, and that Mq is an approximation to the M+ 
solution. This Mq cannot be the real mass of the 
cloud because (i) it is, at best, its upper limit, (ii) it 
is on the unstable branch for the solution of M, and 

(iii) connected to the first reason, it can overesti¬ 
mate the cloud mass (mostly likely M_) by orders 
of magnitude, since Mq/M_ = l/a_ ^ 1 when 

/3< 1. 

Point (4) highlights the necessity of including external 
forces in the analysis of the stability of clouds, no matter 
how small the external forces may be. This point is often 
overlooked in the literature. The root of this ignorance 
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is the intuition that self-gravity is the dominant binding 
force in clouds. 

This intuition is fallacious, as one can see by evalu¬ 
ating the relative importance of the two binding forces 
in the Virial theorem, i.e. comparing the external pres¬ 
sure term and the self-gravity term aGM'^/Rc 

in Equation H. Their ratio, (1 — a)la^ is greater than 
1 for a dynamically stable cloud, i.e. external pressure 
predominates, because we have proven that a must be 
smaller than 1/2 to qualify as a stable solution. If the 
ratio is smaller than 1, i.e. if the self-gravity predom¬ 
inates, the value of a cannot be elsewhere but greater 
than 1/2, and we have proven that such a solution is 
dynamically unstable. 


2.2. Tidal Terms 

A galactic nucleus, which is the focus of this work, 
usually has a SMBH sitting at the centre with a dense 
(nuclear) star cluster surrounding it. Their presence in¬ 
troduces in our problem an additional ingredient of tidal 
fields, which will result in new terms to be added to the 
Virial theorem. 

Besides feeling tidal forces, our test cloud is very likely 
to have a certain degree of rotation as well. This is so 
because its orbital motion in the galactic nucleus leads 
to periodic tidal perturbations, which in the long run 
tends to synchronize the spin o f the cloud with its orbital 
period (|Gladman et al.||1996|). At this point, we say 
that the cloud is "tidally locked”. In this subsection, 
we will include the effects of rotation in our analysis by 
exploring two representative cases: (i) When the cloud is 
non-rotating with respect to the rest frame - the observer 
on the Earth - and (ii) when it is already tidally locked. 

In the following, we will show that (a) the Virial the¬ 
orem without tidal forces, (b) the Virial theorem with 
tidal forces but no rotation, and (c) the Virial theorem 
with tidal forces and a tidally-locked cloud can all be 
merged into a new, universal theorem. 

We start our analysis by adding in Equation Q an 
external gravitational potential $extj accounting foimoth 
the SMBH and the nuclear star cluster. The equation 
of motion of a fluid element inside our test cloud hence 
becomes 


pR =-VP - - pVfhext- (8) 

In the last equation we have introduced R as a position 
vector for the fluid element with its origin at the SMBH, 
which we assume to be at the galactic centre. We now 
introduce a new position vector, D, which has the same 
origin but points at the Centre-of-Mass (CoM) of the 
cloud. The difference of the two, R — D, is a vector that 
points from the CoM to the fluid element of the cloud, 
which we name r. Figure illustrates the definitions of 
the vectors. These new vectors allow us to re-write the 
last equation as follows: 


pf = _VP - + p[-V$ext (R) - D]. (9) 

^-V-' 

tidal term, pat 

In the brackets we can now easily identify - thanks to 
the introduction of D and R - the tidal acceleration a^. 



Fig. 3. — Reference frame for a fluid element inside a non-rotating 
cloud which is orbiting a SMBH. Starting from the SMBH, the 
vector R points to the fluid element and D points to the CoM 
of the cloud. The coordinates are chosen in a way such that the 
origin coincides with the CoM of the cloud, the fc-axis points in the 
same direction as D, and the z-SLX.is is perpendicular to the orbital 
plane of the cloud. In this coordinate system, the vector R — D is 
pointing from the CoM of the cloud to the fluid element, which we 
denote as r, and its three components are (x, y, z). 

We now proceed the same way as we did in Section [24] 
- take the dot product with respect to r and integrate 
Equation ® over the volume of the cloud. The inte¬ 
gration is different depending on whether the cloud is 
rotating or not. 


2 . 2 . 1 . Non-rotating Clouds 

In this case, the orientation of r is fixed with respect 
to an observer standing on the earth, so that the inte¬ 
gration is almost straightforward. As a matter of f act, 
the integration has already been done in Section 2.1 (see 
Equations 1^ to except for the last tidal term p&t. 

To complete tne integration, we need to know the three 
components of the vector a^, so we have to define our 
working coordinates. As is illustrated in Figure [^ we 
choose the origin of coordinates at the CoM of the Uoud. 
The x-axis is chosen to be aligned with D and the z-axis 
aligned with the orbital angular momentum of the cloud. 
Hence, r = (x, y, z). We define ojc to be the angular 
velocity of a circular orbit at a distance of D from the 
SMBH, so that D = 0, 0) if the nuclear star 

cluster is spherically symmetric. 

Civen that |r| |D| in the system of our interest, we 

do a linear expansion for the components of and find 


at = (xT, -yujl, -zujI). 


( 10 ) 


Here we have introduced T = —Ddoj'^/dD, the tidal ac- 
celeration per unit length in the radial direction (e.g. 
Stark & Blitz|1978). The two components with negative 
signs correspond to compressive tides: A fluid element 
moving away from the CoM of the cloud will experience 
a restoring acceleration in the y and z directions. Civen 
that T is positive in galaxies (as we will see later), the 
acceleration in the x direction is pointing away from the 
CoM. 

Therefore, we have seen that the tidal forces acting on 
the cloud is not spherically symmetric. Consequently, a 
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Fig. 4.— Reference frame for a fluid element in a rotating cloud. 
The two position vectors, R and D, and the origm of the coor¬ 
dinates are defined in the same way as in Figure^ but the ori¬ 
entations of the coordinates are not. Here, the tL!-axis is deflned 
by the rotating axis of the cloud, and the other two orthogonal 
axes, namely u- and u-axes, are co-rotating with the cloud. Now 
we denote the vector R — D as f, to distinguish it from the r vector 
in the non-rotating case, and in the current coordinate system f 
has three components (u, u, w). If the cloud is tidally locked, the 
tu-axis will be perpendicular to the orbital plane of the cloud. 


f = 


(ii, V, w) 


In the rotating frame centrifugal 


(13) 


where the acceleration acen comes from the usual con¬ 
cept of a centrifugal force. Let us consider now the con¬ 
sequences for the integration: 

• The first acceleration corresponds to the motion of 
a fluid element with respect to the chosen frame. 
Because we have chosen the co-rotating frame, the 
fluid elements of the cloud are “frozen” with respect 
to the frame, so an integration of the dot product 
of this acceleration with f will vanish, exactly as in 
Equation <§■ 

• As for the centrifugal acceleration, the integral 
Jy f • Siccndm does not vanish because the product 
f • acen is always positive, except for those fluids ly¬ 
ing on the rotation axis of the cloud, in which case 
f • Hcen will be zero. 


Aft er these considerations, we can rearrange Equa¬ 
tion (12), moving all the terms with non-vanishing in- 
tegrals to the RHS, and derive a new, more convenient 
form for the equation of motion: 


cloud in hydrostatic equilibrium is not strictly spherical. 
However, for mathematical simplicity, we inherit spheric¬ 
ity for the later analysis and we deem this approximation 
valid to first-order. This small sacrifice of mathematical 
accuracy will reward us with useful physical insights, as 
we will see. Now using the three components of at and 
assuming sphericity, we derive 


/ v-atdm = bR^M{T -2ujI), (11) 

Jv 

where 6=1/5 for spherical, homogeneous clouds. 

2.2.2. Rotating, Tidally-locked Clouds 

The second representative case in our study is a cloud 
that is tidally locked in its orbit around the SMBH. In 
this case, the integral of the LHS of Equation (1^, i.e. 
Equation (i): is no longer zero. To address the calcula¬ 
tion, we now choose a more convenient frame, the one 
centered on the CoM of the cloud and co-rotating with 
it. Figure [^illustrates this co-rotating frame. 

From the CoM to a fluid element in the cloud, the 
position vector is defined to be f := ue^ + v ey + w e^, 
with Bi (i = u, V, w) the unity vector of components. As 
in the previous non-rotating case, the relation R = D-f-f 
still holds. Therefore, starting from Equation (j^, we 
derive the equation of motion for the considerea fluid 
element as: 


p {it, V, w) = -VP - + P (a* + Seen)- (14) 

We are now left with only one task, namely the integra¬ 
tion of the tidal term Jy f • (a^ -|- acen) dm, because we 
already know the outcome of the other integrals from the 
previous subsection. 

If the cloud is tidally-locked and on a circular orbit 
around the SMBH, the angular velocity (or spin) vector 
= (0, 0, oje) will be constant, and the centrifugal ac¬ 
celeration can be calculated with acen = x (FJ x f) = 
(uw/, uw/, 0). After some algebra we find 


f f ■ {at + accn)dm = bR^MT, 

JV 


(15) 


where, again, we assume sphericity for the cloud. 

2.3. The Extended Virial Theorem 

Before we continue, we recapitulate what we have got 
so far: (i) The commonly-used virial mass Mq is in¬ 
evitably an overestimation of the mass of a dynamically 
stable cloud, (ii) The tidal forces lead to an extra term 
in the derivation of the Virial theorem, and we shall 
call this new theorem the extended Virial theorem (EVT 
hereafter), (iii) We evaluated two representative cases to 
calculate this extra ter m, n amely a non-rotating cloud, 
which led to Equation (11), and a tidally-locked cloud. 


which resulted in Equation (15). 

It is easy to see that the EVT takes the form: 


pf =-VP - pV$c + pat. (12) 

This equation is similar to Equation (j^, except that 
the unit vectors that we choose for f are rotating respect 
to an observer on Earth. As a result of this choice, f is 
a combination of two accelerations, i.e. 


47rP/Ps = - aGM^/Ry, (16) 

if we define a', the effective velocity dispersion, in the 
following way: 

ct' 2 =cr^-H6P^(r-2a;/)/3 (17) 
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for non-rotating clouds and 


= a^ + bR^T/3 


(18) 


for tidally-locked ones. Now it is clear that an exter¬ 
nal tidal field effectively changes the internal energy of a 
cloud, in other words, it changes the ability of the cloud 
to expand or to contract. _ 

The discussion presented in Section |2.1| must be re¬ 
evaluated in the light of our newly derived Virial the¬ 
orem. Here we highlight two cru cial points, (i) It can 
readily be seen that Equation (16 1 is structurally identi¬ 
cal to Equation (§. The only ditterence is the newly de¬ 
fined cr', which equivalently introduces a new virial mass, 
Mq = 2>R(j'^ / [aG). Correspondingly, the dimensionless 
mass should be redefined as a' = M/Mg and the dimen¬ 
sionless surface pressure as f3' = APsV /so that 
the EVT retains a dimensionless form of /3' = 4q;'( 1 —o'), 
(ii) The ETV is a quadratic function of M, same as 
the conventional Virial theorem. Therefore, we can 
prove that given f3', the dynamically-stable solution is 
M'_ = Qf'_Mg, where a(_ < (1/2). Again, this solution 
stands for a cloud bounded mainly by the external forces. 

We now rewrite a' in a different way, so that it is easier 
to understand what the implications of our EVT are. If 
the central SMBH has a mass M, and the nuc lear star 
cluster follow s a density profile p* oc D~~> (see Amaro- 
Seoane 2012 for a review), the angular velocity Wc of a 
circular orbit can be written as 


2 G[M. 
uj: = - 


M*(D)] 


£>3 


(19) 


where M^,{D) oc is the stellar mass enclosed by the 

orbit. It follows that 


..r'2 _ 

a = a 


bR^ 


.black hole 

for the non-rotating clouds 


GM, 


- (2 - 7 ) 


GM^{D) 


stellar component 


( 20 ) 


^'2 2 
(7 — (7 


bR^ 


3GM, GM^{D) 


D3 


■7" 


£)3 


( 21 ) 


for tidally-locked clouds. It is clear that, as marked in 
Equation (20), two components contribute to the tidal 
effects. The first contribution is directly linked to the 
mass of the SMBH, and the second one to the mass and 
shape of the stellar system. 

A close inspection of the last two equations reveals two 
pieces of new information: 

(i) If the cloud is non-rotating, then we could have a' < 
tr, i.e. when the SMBH term of Equation (20) becomes 


smaller than the stellar one. In this case the tidal force is 
effectively reducing the internal energy of the cloud. This 
effect is coming from the stellar component, especially its 
tidal forces in the tangential direction which are acting in 
a compressive way on to the cloud. On the other hand, 
we can have the inverse situation, in which a < a', i.e. 
when the cloud is so close to the SMBH so that the role 


of the stellar system can be neglected (very small D). 
In this latter situation, the tides, induced mainly by the 
SMBH, essentially are tearing the cloud apart, increasing 
the internal energy of the cloud. 

(ii) In the tidally-locked case, we have a similar situ- 
ation, becaus e 7 is not necessarily always positive (e.g. 


Merritt 2010). Contrasting the last two equations, one 
can see that a cloud that is tidally-locked always has a 
larger a' than that of a non-rotating one. This excess of 
effective internal energy comes from the additional cen¬ 
trifugal acceleration. 


3. APPLICATIONS 

In this section we apply the stable solution from the 
EVT, i.e. M'_ = a'_MQ, to infer the physical properties 
of observed gas clumps in two types of galactic nuclei, 
namely AGNs and our GC. As we will see, using M'_ in¬ 
stead of the conventional virial mass Mq resolves a series 
of controversies and provides interesting implications for 
the origin and evolution of these clumps. 

From now on we will call an observed object a “clump”, 
to distinguish it from the idealized, theoretical concept 
of a “cloud”. This distinction is meaningful because it 
draws attention to the limitations of the observational 
data - while clouds have well-defined shapes and physical 
parameters, clumps do not have clearcut boundaries, nor 
are all their physical quantities observable. 

In each of the three subsections that follow we (i) first 
introduce the current controversy linked to the usage of 
the conventional virial mass and (ii) resolve the contro¬ 
versy by applying our EVT. 


3.1. Clumps in AGN Tori 

Controversy: Analysis of the X-ray variability of 

AGNs indicates the existence of a population of dusty 



clump. We summarize here the typical values for them: 
(i) A size of Rc — 2.0 x cm, (ii) a Hydrogen (num¬ 
ber) density of about uh — 3.5 x 10® cm“3^ (iii) a mass 
for the central SMBH of M, ~ 4.2 x 10^ Mq, (iv) a dis¬ 
tance of 13 ~ 5.0 X 10^^ cm from the SMBH , and (v) 
a bolometric luminosity of L ~ 1.6 x 10 "^'* erg s“^ for 
the AGN. From the size and density, we derive a typical 
mass of M = 10“® Mq for a spherical and homogeneous 
clump. 

Naively, one would expect a clump to be tidally dis¬ 
solved if the tidal force from the SMBH acting on to the 
clump exceeds its self-gravity. Hence, it would seem that 
the density of a clumps should be at least 


ntide - M,/(mpI3®). 


( 22 ) 


From the quoted typical values and a proton mass of 
nip = 1.67 X 10“^"* g, this critical density is ntide — 
4.0 X 10^^ cm“3^ a thousand ti mes above nn. And yet the 
clumps exist. This paradox led Markowitz et al. (2014) to 


discuss possible stabilizing mechanisms, including mag¬ 
netic field and external pressure. 
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Solution: The aforementioned clump - whose self¬ 

gravity is weaker than the external tidal forces but re¬ 
mains bound and stable - is exactly the type of cloud 
we have found as the stable solution of our EVT (Sec¬ 
tion 2.3). Our EVT predicts that the clump has a non- 
negligible surface pressure P 5 , and its mass is q;'_Mq 
where a'_ < (1/2). Therefore, we equate a'_MQ to 
10 “® Mq, the typical mass of the clump, and we pro¬ 
ceed to calculate a'_ and then Pg. 

To do so, we have all the data we need except those 
of (T, 7 , and M^{D). However, it is relatively str aight¬ 
forward to obtain a well-established value for a. IKro-l 
lik & Lepp (1989) proved that the typical gas temper¬ 
ature in a rnolecular cloud in an AGN torus is about 
10^ K. This corresponds to a typical internal velocity of 
tr « 3 km s“^. 

As for 7 and M*(P), which determine the value of cr', 
we notice that M^{D) M, at the distance D ~ 0.2 
pc where the clumps are observed. Therefore, we can 
neglect the stellar components, i.e. those determined by 
M*(P)) and 7 , in the analysis of the tidal terms. 

Although the SMBH term predominates in the brack¬ 
ets of Equations (20) and (21), the two terms in the 
brackets are multiplied by Pj/This multiplication leads 
to a negligibly small contribution compared to cr^, so es¬ 
sentially tj' « cr « 3 km s“^. It is important to know 
in advance that (i) the tidal term from the SMB H is not 
always negligible, as we will see in Section 3.2 and (ii) 
neither is the ste llar component always negligible, as is 
shown in Section 13.31 

We have specified the values of all the terms entering 
the EVT, so we can solve a'_ and Pg. Erom a' and Rc 
we derive M/ ~ 0.67 Mq. It follows that a'_ = M/Mq ~ 
1.5 X 10“®. Since the value we have just derived for a'_ is 
smaller than 1/2, the mass M = 10“^ Mq qualifies as a 
stable solution in the context of Equation (16). Since we 
know that a'_ must satisfy the relationship)^ = 4 q!'_ (1 — 
a'-), we derive that /3' ~ 4a'_ ~ 6.0 x 10“®, and from 
the definition of /3' (Section 2.3), we finally find Pg ~ 
5.4 X 10“^ erg cm“^. 

Now that we have the value for the surface pressure 
of a stable cloud in our theory, we calculate with ob¬ 
served data the typical external pressure of the environ¬ 
ment where real clumps are found and compare the two. 

The main source of external pressure is the turbulent 
gas surrounding the clump. Estimating its magnitude 
requires knowledge of the mean density next and the tur¬ 
bulent velocity Ve^t of such external gas medium. To 
estimate Uext, we recall the f ollowing two facts, (i) AGN 
tori are dusty (Netzer 2015) and dust grains will be de¬ 
stroyed if the colhsio nal velocity induced by turbulence 
exceeds 120 km s“^ ( Shull|[T978 1 . (ii) The turbulent ve¬ 
locity must also be comparable to the Keplerian velocity, 
about 10^ to 10^ km s“^ depending on M, and D, to 
maintain t he empirically required ge ometric thickness of 
the torus ( Krolik fc Begelman||1988 ). Taking both facts 
into account, we adopt Vext = 100 km s“^. 

Eor next, we know that the integrated Hydrogen col¬ 
umn density in the mid-plane of the torus is an observ¬ 
able, and the ty pical value is Nff ~ lO^'^cm”^ (Buch¬ 


ner et al. 2015). If we estimate the width AD of the 


one of the torus - we can derive next approximately by 
NH/{AD) . To estimate AP, we use two pieces of empir¬ 
ical information: (i) there is a correlation between AGN 
lumin osity and the radius of the inner boundary of the 
torus (|Nenkova et al.|2008|), according to which the inner 
boundary is at 0.2 pc tor our typical AGN with a lumi¬ 
nosity of P ~ 1 . 6 x 10 "^'* erg s“^, and (ii) the outer bound- 
ary is typically 10 times more distant than the inner one 
(Elitzur [M07 I, so it is at 2 pc in our case. As a result. 


we find AD ~ 1.8 pc, and hence next — 1.8 x 10® cm“®. 

With the numbers that we have derived for z;ext and 
next, we proceed to calculate the external pressure with 
Text — ^ext , and the result is 2.9X 10 ® erg cm 
We find that Pgxt is consistent with Pg within a fac¬ 
tor two. This agreement supports our proposal that the 
clumps are confined by an external pressure and are dy¬ 
namically stable. 

We note that the radiation from the central AGN can¬ 
not be another source of external pressure. This is so, 
because a dusty cloud exposed to AGN irradiation will 


evaporate within one orbital period (Pier h Voit 1995, 
Namekata et al.|[2014 1. The radiative pressure would be 


Prad = LI{At:D‘^c), (23) 

with c the speed of light. We can see that this value is 
5.1 X 10“® erg cm“®, more than two orders of magnitude 
higher than the surface pressure Pg derived above for the 
clump. Hence, for the clump to not be dissolved, it must 
be shielded during most of its orbital phase by the gas or 
other clumps loca ted closer to the central SMBH, as has 
been sugges ted by Krolik & Begelman (1988); Namekata 
et al. I (120141). Our result, hence, supports this shielding. 


3.2. Clumps in the CND of the Milky Way 

Controversy: Molecular clumps have been detected 
in the centers of quiescent galaxies as well (Hunt fc 


Malkan 


2004| |Mazzalay et al. 20131 |Davies et a,l.||2U14|). 
In the Go, in particular, tens ot molecular cluinps re¬ 
side in the GND, a ring-like (narrow projected width) 
structure surrounding Sgr A*. The observable parame¬ 
ters and their typical values a re D = 1.8 pc. = 0.2 5 
pc, and (T = 11 km s“^ (from Christopher et al. 2005) 


and the central S MBH has a well-es tablished iiiass of 
M, « 4 X 10® Mq (iGenzel et al.||2010|). 

There exists two ditterent methods to derive the masses 
of these clumps. However, both lead to dilemmas that 
seem to invalidate the method. We now explain the 
methods and the problems that emerge with them. 

(i) The typical mass of the clump according to the 


3.4 X 10^ Mo falso see Shukla et al.||2004 

Christo- 

pher et al. 1120051 Montero-Gastaho et al.||20U9 

ITsuboi 

k Miyazakil 2012 

). A clump with this mass has a den- 


torus - measured from the inner boundary to the outer 


, comparable to the tidal density 
ritide — 2.8 X 10^ cm“® at the observed distance {D), so 
the clump can withstand the tidal force. This argument 
has been used to support the validity of taking Mq as the 
real clump mass. However, this mass has two problems, 
(a) In the context of the conventional Virial theorem, 
such clumps are gravitationally un stabl e, as discussed in 
points 3 and 4 of the list in Section [2T] They will disap¬ 
pear on a free-fall timescale (about hP yr), shorter than 










































































the required time to form an axisymmetric CND (about 
10^ yr, comparable to the orbital period). In the EVT, 
however, tidally-locked clumps with this mass are stable, 
as we will see later, (b) Nonetheless, Mq is likely to be 
ruled out, because if this were indeed the real mass of the 
clumps, then the CND would have a total mass of about 
10® Mq. Such a heavy disk would inevitably destroy the 
observed young st ellar disk within 0.5 pc around Sgr A* 


Subr et al. 

2009 

I. 

(ii) Anot 

ler way ot deriving the mass of clumps 


is based on the assumption of local thermal equilib¬ 
rium (LTE, and in this case we call the mass derived 
Mute)- Several authors obtained a similar result, typ¬ 
ically Mlte — 2.6 X 10^ Mp), one order of magn itude 


smaller than before (from Tsuboi & Miyazaki |2012| also 

see IMarr et al. 119931 IMarshall et al.| 

19951 |(Jhristoplier| 

et al. 2005| Liu et al.||2013| Requena-l 

forres et al.||2012). 

in this case, Mlte does not lead to t 

he aforementioned 


problem related to the existence of the young stellar disk 
in our GC. However, according to the naive tidal limit 
criterion, these clumps will disappear within a relatively 
short time, in any case shorter than the required time 
to form an axisymmetric CND, as has been noticed by 
these authors. 

Solution: We will see that in the context of the EVT, 
Mq, as derived in (i), is unstable for non-rotating clumps 
but stable for tidally-locked ones. Nevertheless this does 
not change the fact that this mass leads to a scenario 
ruled out by observations. On the other hand, a clump 
with mass Mlte, as derived in (ii) can be stable. Since 
this mass does not lead to contradictions with observa¬ 
tions of the GC, it is a realistic value for clumps in the 
CND. 

We now derive the value of Ps as provided by the EVT, 
and then we will compare it with the environmental pres¬ 
sures independently derived in other means, to prove that 
the clump with a mass of Mlte is in hydrostatic equilib¬ 
rium. 

The GC nuclear cluster has 7 = 1.75, and the enclosed 
stellar mass is approximately 


M^D) ~ 1.6 X (24) 


T his stellar distribution is adopted from the early w ork 
of Vollmer & Duschl (2000); Christopher et al.|(2005) for 
the sake of comparison, but it is not signihcantly ditterent 
from t he distribution der ived from more recent observa¬ 
tions ( Genzel et al.|[2M0 ). With this stellar component, 
we find a' = 14 km s“^ for a non-rotating clump and 
tr' = 19 km s“^ for a tidally-locked one. 

Eor the non-rotating clump, we hence derive Mg ~ 
5.9 X 10^ Mq, = Mlte/M^ 0.044 « 1, /3' ~ 
4a'_ ~ 0.18, and finally Ps — 6.7 x 10“® erg cm“®. Eor 


the tidally-locked one, Mg = 1.1 x 10® Mq, a'_ ~ 0.024, 
/?' ~ 0.094, and Ps — 1-2 x 10“® erg cm“®. In the two 
cases, we have a'_ < (1/2), so we can conclude that M = 
Mlte is a stable solution. We also note that M = Mg is 
an unstable solution for non-rotating clumps, but it could 
be a stable solution for tidally-locked ones, although, as 
mentioned before, this mass is excluded because of other 
arguments not related to the EVT. 

Balancing this surface pressure requires an environ¬ 


mental pressure of at least the same order of magnitude. 
In the GG we have three possibilities: First, the hot. 


of 10^ cm ® (Jackson et al 

19931 and a temperature of 

7,000 K (Kobe^ & Goss 

99^3), which provides a ther- 


weak to explain the magnitude of Ps- Any derivation of 
the mass of the clumps based on this pressure will lead to 


a too small value , such as the 15 Mq derived by (|Vollmer 
fc Duschl|2001a I. Second, at a distance of D = 1 pc, the 


ram pressure caused by the stellar winds from t he ob¬ 
served young stars is a bout 1.1 x 10“^ erg cm“® (lYusef- 

. TM 


Zadeh fc Wardle|1993 ), which is still too weak. The last 


possibility otters a plausible solution: The molecular gas 


lent velocity of 55 km s 


-1 


(Giisten et al. 

19871) and 

cm“'^ (Genze 

1 et al.||1985 


5 X 10“® erg cm“®, consistent with Ps as derived above. 
Hence, the clumps in the GND can be stabilized by the 
inter-clump turbulent gas. 

The turbulence of the inter-clump gas has been postu¬ 
lated to be generated by a recent outburst of Sgr A* and 
continuously replenished b y the dissipation of the differ¬ 
ential rotation of the GND (Giisten et al.|19^). We have 


now a number of evidences suggesting that the outburst 


took place about 6 Myr ago and th( 
ity was as high as 10^® erg s“^ (e.g. 

5 bolometric luminos- 

Nayakshin & Cuadra 

20051 ISu et al.|2010l|Bland-Hawthorn et al. 120131 

Amaro- 

Seoane & Ghen|2014 Chen & Amaro-Seoane|201'^ 

ll|2015p. 


on the inner rim of the GND (at about D = 1 pc), of 
about 3 X 10“® erg cm“® (see Equation ). This value 
is comparable to Ps as well as to the turWlent gas pres¬ 
sure. Therefore, we infer that during the outburst of Sgr 
A*, the inner rim of the GND was pressurized by the ra¬ 
diation and that the pressure propagated into the CND 
to stabilize the clumps in the context of the EVT. 


3.3. Clumps in the CMZ 

Cont rove rsy: The methods (i) and (ii) described in 
Section [3?^ have been applied to the clumps in the CMZ 
of the Milky Way to derive their masses. The result is 
that the vir ial mass Mg is on average 10 times large r 
than Mlte ( Miyazaki fc Tsub^|2000 Oka et aL]|200i ), 
a difference as large a s we h ave seen m Section |3.2[ 


Miyazaki & Tsuboif (2000) investigated the possibility 
that the clumps m the GMZ are confined by the hot, 
inter-clump ISM. However, the pressure of the hot ISM, 
obtained from the observations of the 6.7 keV iron lines, 
is ten times smaller than the surface pressure that one 
derives based on the conventional Virial theorem (not the 
EVT). Therefore, the external pressure could not explain 
the discrepancy. 

Partial solution: We adopt now D — 100 pc, Rc = 10 
pc, and a = 10 km s“^ as the fiducial values (from Oka 
et al. ]20M| ), out we note that the CMZ spans a large ra- 
dial range, and the clumps in it have a broad distribution 
in size as well as in internal velocity dispersion. 

In Section |2.3| we addressed two particular cases, 
namely tidally-locked and non-rotating clumps. In the 
following, we will show that it is the case of non-rotating 
clumps what (partially) solves the present problem re- 
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lated to the clumps in the CMZ. 

It does so, because the effective velocity dispersion a' 
is smaller than the conventional one a. Consequently, 
applying the Virial theorem without accounting for the 
tidal forces will lead to an overestimation of the surface 
pressure of the clump. We can see this by comparing 
the conventional an d the extended Virial theorems, i.e. 
Equations [6| and [TG] Physically, this can be envisaged as 
the result iff ignoring the contribution of the compressive 
tides in the stabilization o f th e clumps. 

According to Equations ([20l) and given that 2—7 > 0 in 
the Milky Way ( Genzel et al.j2010 ), the situation o' < a 
will occur, and it happens when the (effective) stellar 
mass (2 — 7)M*(iD) enclosed by the orbit of the clump 
exceeds the mass M, of the SMBH. The corresponding 
region is at I? > 6.3 pc, which encompasses the entire 
CMZ. 


To be more quantitative, we scale Equation (20) using 
the parameters releva nt to the clumps in the CMZ (also 
from Oka et ar||2001 1 and we find 


~ cr^ - (4.3 km s-i)2 b 


Rr 


10 pc J \100 pc 


D 


-1.75 


(25) 

Given u ~ 10 km s“^ and D ^ 10^ pc, the last equation 
suggests that if a clump is bigger than 10 pc, its surface 
pressure will be significantly weaker than what was de¬ 
rived in the early work based on the conventional Virial 

theorem. _ 

In a survey of the CMZ ( Oka et al.|2001 ), 67 out of 165 
molecular clouds have been observed to have a size larger 
than 10 pc. Limited to these observations, our EVT par¬ 
tially resolves the controversy. Other explanation must 
be sought to fully address the problem (see Section]^ for 
discussion). 

4. DISCUSSIONS AND CONCLUSION 


homogeneity in the ISM can either compress or decom¬ 
press a test cloud. The fastest relative velocity between 
the cloud and the ISM is the Keplerian velocity, there¬ 
fore the shortest timescale for the cloud to return to the 
inhomogeneous region is the orbital period, torb- 
After one perturbation, the cloud requires at least the 
sound crossing timescale tsc — 2i?c/<7 to relax and re¬ 
store equilibrium. So our assumption of hydrostatic equi¬ 
librium is valid if these two criteria are simultaneously 
met: (i) Gc < AoU and (ii) tsc < torb- We now use these 
two criteria to assess the destabilizing effects in the three 
different media that we have considered in this paper. 

(a) Eor a cloud in AGN torus, since we know the typical 
values for cr, Rc, D. M,. we find that Af ::::: 42 yr, torb — 
6.3 X 10^ yr. Since Krolik & Begelman (1988) proved that 


tcoii — torb, we conclude that hydrostatic equilibrium can 
be achieved in AGN tori. 

(b) For a cloud in the CND, we find tgc — 4.5 x 10^ yr 
and torb — 1.1 X 10^ yr using the observational data. 
From the numerical simulations of IVollmer fc DuscEH 
(2002), we have that tcoii — 2 Myr. Hence, hydrostatic 
equilibrium can be achieved in the CND. 

(c) For a cloud in the CMZ, we derive tgc — 2 Myr and 
torb — 5.9 Myr according to the observational data. For 
tcoii) we n otice that the clu mps cover about 50% of the 
disk area (Oka et ar]|2001|), so we use the formula for 


tcoii derived m Bally et al. 1|1988 1, which is a function of 


the covering factor, and we find icoU — 4 Myr for a cloud 
at D = 100 pc. In this region, the three timescales (i.e. 
tsc, Aoii, and torb) are of the same order of magnitude. 
As a result, the clouds may not have had enough time 
to fully restore equilibrium. This could be the reason 
why in Section |3.3| we only could explain a fraction of 
the clumps in the CMZ. Indeed, there is evidence that 
at the inner rim of the CMZ, where both tcoii and torb 
are relatively short, the molecular clouds may have re¬ 
cently expe rienced collisions (e.g. Rodriguez-Fernandez 
et al.| |20(i6 |. 


In this paper we analyze the stability of the gas clouds 
in the central 1 —10^ pc of a galaxy, in the region around 
the SMBH and the nuclear star cluster. We have s hown 
that the external forces - external pressure (Section |2.1[ ) 
and tidal forces (Section |2.2[ ) - play a crucial role in con¬ 
fining and stabilizing the clouds. Based on this under¬ 
standing, we formulated an extended Virial theorem (or 
EVT) to identify the correct ranges of p hysical quan¬ 
tities that lead to stability (Section 2.3). We applied 
our EVT to model observational data and have solved 
practical problems related to the sta bility of those gas 
clumps detected in A GN tori (Section [3.l[ ), the CND in 
the GC (Sec tion 3.2), and the CMZ of the Milky Way 


(Section 

problem). 


in this last case we partially resolved the 


4.1. Validity of the Assumption of Static Equilibrium 

In the previous sections we have relied on the assump¬ 
tion that the clouds are in hydrostatic equilibrium. How¬ 
ever, there are two kinds of perturbations that poten¬ 
tially could invalidate this assumption, (i) Collisions 
with other clouds could drive a test cloud out of equi¬ 
librium, if the cloud number density is high enough. We 
define tcoU as the collisional timescale - the mean times- 
pan during consecutive collisions, (ii) The effect of in¬ 


4.2. Conclusion 

We have proven that the observed gas clumps in AGNs 
tori and in the CND of the Milky Way (and partially 
in the CMZ) are not self-gravitating, but instead, they 
are stabilized by external pressure and tidal forces. The 
standard Jeans instability produces only self-gravitating 
clouds, and therefore it cannot account for the origin of 
these observed clumps. 

The noticeably large internal energy relative to the self¬ 
gravity, as we found for these clumps, is more consistent 
with a non-conventional formation mechanism, such as 


the collisional fragmentation and agglo meration (Krolik 
& Begelman 1988 Vollmer et al. 2004). This discovery 
supports the picture that clump-clump collision drives 
gas inflow to the central pc to power AGNs . 

In the future, application of the EVT to individual 
clumps (i.e. observed clouds) potentially can provide us 
with a mapping of clump masses and environmental pres¬ 
sure around supermassive black holes in galactic nuclei. 
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